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ABSTRACT 

We present new 230 GHz Submillimeter Array observations of the candidate first hydrostatic core 
Per-Bolo 58. We report the detection of a 1.3 mm continuum source and a bipolar molecular outflow, 
both centered on the position of the candidate first hydrostatic core. The continuum detection has 
a total flux density of 26.6 ± 4.0 mJy, from which we calculate a total (gas and dust) mass of 
0.11 ± 0.05 M Q and a mean number density of 2.0 ± 1.6 x 10 7 cm -3 . There is some evidence for the 
existence of an unresolved component in the continuum detection, but longer-baseline observations 
are required in order to confirm the presence of this component and determine whether its origin lies 
in a circumstellar disk or in the dense inner envelope. The bipolar molecular outflow is observed along 
a nearly due east-west axis. The outflow is slow (characteristic velocity of 2.9 km s _1 ), shows a jet- 
like morphology (opening semi-angles for both lobes), and extends to the edges of the primary 
beam. We calculate the kinematic and dynamic properties of the outflow in the standard manner 
and compare them to several other protostars and candidate first hydrostatic cores with similarly low 
luminosities. We discuss the evidence both in support of and against the possibility that Per-Bolo 
58 is a first hydrostatic core, and we outline future work needed to further evaluate the evolutionary 
status of this object. 

Subject headings: stars: formation - stars: low-mass - ISM: individual objects (Per-Bolo 58) - ISM: 
jets and outflows 



1. INTRODUCTION 

Low-mass stars form from the gravitational collapse of 
dense molecular cloud cores (e.g., Shu, Adams, & Lizano 
1987; McKee & Ostriker 2007). Although many open 
questions remain, the evolutionary stages in this process 
are generally well understood. Nearly all of the main 
stages have been observed in large samples and compared 
to theoretical predictions, including starless (prestellar) 
cores, embedded protostars, T Tauri stars, and pre-main 
sequence stars. However, one stage has still not yet been 
definitively observed: the first hydrostatic core (here- 
after FHSC). First predicted by Larson (1969), it ex- 
ists between the starless and protostellar phases. As the 
molecular cloud core collapses, a central hydrostatic ob- 
ject (the FHSC) forms once the central density increases 
to the point where the central region becomes opaque 
to radiation (p c > 10~ 13 g cm~ 3 ; Larson 1969), render- 
ing the collapse adiabatic rather than isothermal. This 
object continues to accrete from the surrounding core 
and both its mass and central temperature increase with 
time. Once the temperature reaches ~ 2000 K the energy 
liberated by accreting material dissociates H2, prevent- 
ing the temperature from continuing to rise to sufficiently 
balance gravity. At this point the second collapse is ini- 
tiated, eventually leading to the formation of the second 
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hydrostatic core, more commonly referred to as the pro- 
tostar. 

Numerous theoretical studies have investigated the 
physical properties of first hydrostatic cores. The max- 
imum mass before the onset of the second collapse is 
generally - 0.04 - 0.05 M (Boss & Yorke 1995; Ma- 
sunaga, Miyama, & Inutsuka 1998; Omukai 2007; To- 
mida et al. 2010), although a larger range (0.01 — 0.1 
M Q ) is found when the effects of rotation are included 
(Saigo & Tomisaka 2006; Saigo, Tomisaka, & Matsumoto 
2008). The lifetime of the FHSC is estimated to be short 
compared to the embedded phase duration of 0.54 Myr 
(Evans et al. 2009), but specific studies find a large range 
of possible lifetimes between ~ 5 x 10 2 — 5 x 10 4 yr (Boss 
& Yorke 1995; Omukai 2007; Saigo, Tomisaka, & Mat- 
sumoto 2008; Tomida et al. 2010). Most of the varia- 
tion is due to different models and assumptions concern- 
ing the mass accretion rate from the surrounding core. 
The predicted internal luminosity 6 of the FHSC is also 
very low but covers a large range (10 -4 — 10 _1 L Q ; Ma- 
sunaga, Miyama, & Inutsuka 1998; Omukai 2007), with 
most of the variation again due to different mass accre- 
tion rates from the surrounding core leading to differ- 
ent accretion luminosities. The radius is typically ~ 5 
AU (Masaunaga, Miyama, & Inutsuka 1998), but can be 
much larger (> 10 — 20 AU) and exhibit a flattened, disk- 
like morphology when the effects of rotation are included 
(Saigo & Tomisaka 2006; Saigo, Tomisaka, & Matsumoto 
2008). 

The observational characteristics of FHSCs have also 
been predicted by several authors. The emergent spec- 

6 The internal luminosity is the luminosity of the central source 
and excludes luminosity arising from external heating. 
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tral energy distributions (SEDs) of FHSCs embedded in 
their parent cores have been studied by Boss & Yorke 
(1995), Masunaga, Miyama, & Inutsuka (1998), Omukai 
(2007), and Saigo & Tomisaka (2011). While the de- 
tails vary from one study to the next, they all find that 
the radiation emitted by the FHSC is essentially com- 
pletely reprocessed by the dust in the surrounding core 
to 50 — 200 /xm, with the SED characterized by emis- 
sion from 10 — 30 K dust and no observable emission 
below 30 — 50 /im. Omukai (2007) also showed that the 
H2O lines emitted in the accretion shock fall below the 
sensitivities of current space missions but could be de- 
tected in future planned missions. Machida, Inutsuka, & 
Matsumoto (2008) showed that both the first and second 
cores drive outflows. The outflow driven by the first core 
features slow velocities (< 5 km s _1 ) and wide opening- 
angles; in contrast, the outflow driven by the second core 
is faster (~ 30 km s _1 ) and well-collimated. 

On the observational front, several candidate FHSCs 
have been found, although none have been definitively 
identified as true first cores. Belloche et al. (2006) found 
evidence for significant chemical evolution in the dense 
core Chamaeleon-MMSl, identified an associated faint 
infrared source in Spitzer Space Telescope (Werner et 
al. 2004) observations, and noted this core harbors either 
a FHSC or very young protostar. Chen et al. (2010) de- 
tected an outflow driven by the previously unknown core 
L1448 IRS2E in Perseus and showed there is no associ- 
ated Spitzer infrared source down to a luminosity limit of 
~ 0.1 Lq. They suggested this is a strong FHSC candi- 
date, although the outflow is much faster (~ 25 km s _1 ) 
than expected (< 5 km s _1 ; Machida, Inutsuka, & Mat- 
sumoto 2008). X. Chen et al. (2011, in preparation) and 
Pineda et al. (2011) have detected outflows from two ad- 
ditional cores lacking Spitzer infrared detections (CB 17 
MMS and L1451-mm, respectively) and have suggested 
both as first core candidates. Enoch et al. (2010) de- 
tected the 70 /xm signature expected for a FHSC in deep 
Spitzer observations of the dense core Per-Bolo 58 in 
Perseus. However, they also found a corresponding faint 
Spitzer 24 /im detection in existing Spitzer c2d legacy 
data (Evans et al. 2003; 2009) and argued this source 
could either be a FHSC or very low luminosity Class 
protostar. 

Determining the true evolutionary status of these ob- 
jects is of central importance in studies of the earliest 
phases of low-mass star formation. Identifying one or 
more as confirmed FHSCs would verify a theoretical pre- 
diction first made over 40 years ago (Larson 1969), and 
would open the physical properties and characteristics of 
FHSCs to observational study. Furthermore, identifying 
one or more as confirmed very low luminosity protostars 
would have important implications for our current under- 
standing of protostellar luminosities and mass accretion. 
The Spiter c2d legacy program has detected a number 
of very low luminosity objects (VeLLOs), protostars em- 
bedded in dense cores with internal luminosities less than 
or equal to 0.1 L Q but generally above ~ 0.05 L Q (Di 
Francesco et al. 2007; Dunham et al. 2008), many in cores 
previously believed to be starless. Recent results from 
this program have also shown that the protostellar lumi- 
nosity distribution is strongly peaked at low luminosities 
(Enoch et al. 2009a; Evans et al. 2009), confirming and 



aggravating the classic "luminosity problem," whereby 
protostars are underluminous compared to the accretion 
luminosity expected both from theoretical collapse calcu- 
lations and arguments based on the minimum accretion 
rate necessary to form a star within the embedded phase 
duration (Kenyon et al. 1990; 1994; Kenyon & Hartmann 
1995). Confirmed detections of protostars with luminosi- 
ties at or below 0.01 L Q would only further exacerbate 
the luminosity problem. 

In this paper we present new 230 GHz Sub-millimeter 
Array (SMA; Ho et al. 2004) observations of Per-Bolo 
58 and report the discovery of an outflow driven by this 
source. We provide a brief introduction to Per-Bolo 58 
in Sj2] and a description of our observations in Sj3j We 
present our basic results in $4l including the detection of 
the continuum in 8 34. H and the detection of 12 CO J = 2—1 
and its isotopologues in £14.21 We analyze the properties 
of the molecular outflow in <j5l including its morphology 
f £|5.1[) and kinematics f £|5.2[) . Finally, we discuss the im- 
plications of our results in Sj6]and summarize our findings 
in 33 

2. PERSEUS BOLO 58 

Perseus Bolo 58 (R.A. = 03 29 25.7, Deck = +31 
28 16.3; hereafter Per-Bolo 58) is a dense molecular 
cloud core located on the northern edge of the young 
(proto)stellar cluster NGC 1333 in the Perseus Molecu- 
lar Cloud, at an assumed distance of 250 ± 50 pc (Enoch 
et al. 2006), consistent with the VLBI maser parallax 
distance of 235 ± 18 pc for NGC 1333 recently deter- 
mined by Hirota et al. (2008). Per-Bolo 58 was identi- 
fied independently by Hatchell et al. (2005) and Enoch 
et al. (2006) in (sub)millimeter dust continuum emission 
surveys of Perseus using James Clerk Maxwell Telescope 
(JCMT) Submillimeter Common User Bolometer Array 
(SCUBA) 450 and 850 /im observations and Caltech Sub- 
millimeter Observatory (CSO) Bolocam 1.1 mm observa- 
tions, respectively. The total core mass is determined to 
be 0.8 - 1.2 M by Enoch et al. (2006; 2010) and 2.6 
M Q by Hatchell et al. (2007a). Schnee et al. (2010) re- 
ported a detection of Per-Bolo 58 in Combined Array 
for Research in Millimeter Astronomy (CARMA) 3 mm 
continuum observations, and from the detection calcu- 
lated a total core mass of 2.4 M Q . The differences in 
these mass estimates are within the uncertainties intro- 
duced by different gas-to-dust ratio, dust opacity, and 
dust temperature assumptions. 

Per-Bolo 58 was identified as starless by both Enoch et 
al. (2008) and Hatchell et al. (2007a) based on the lack of 
an associated Spitzer infrared source. It was also identi- 
fied as starless by Hatchell et al. (2007b) and Hatchell & 
Dunham (2009) based on the nondetection of an outflow 
in single-dish JCMT 12 CO J = 3 - 2 data. Curtis et 
al. (2010) also reported the nondetection of an outflow 
in their independent JCMT 12 CO J = 3 - 2 data. Arce 
et al. (2010) detected a redshifted COMPLETE 7 Perseus 
Outflow Candidate (CPOC 19) in their single-dish Five 
College Radio Astronomy Observatory (FCRAO) 12 CO 
J = 1 — data beginning ~1' northeast of Per-Bolo 58 
and extending several arcminutes farther to the north- 
east. However, they identified IRAS 03262+3123 as the 

7 Coordinated Molecular Probe Line Extinction Thermal Emis- 
sion; Ridge et al. (2006) 
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TABLE 1 

SMA Observations of Perseus Bolo 58 







V 


Beam FWHM 


Beam PA a 


Bandwidth 


<5V b 


1<t rms c 


Line 




(GHz) 


(Arcseconds) 


(Degrees) 


(GHz) 


(km s" 1 ) 


(mjy beam ) 


12 CO J = 2 - 


1 


230.53797 


2.92 X 2.49 


-86.0 


0.082 


0.53 


42 


13 CO J = 2- 


1 


220.39868 


3.15 x 2.63 


-74.5 


0.082 


0.53 


33 


C ls O J = 2 - 


1 


219.56037 


3.16 x 2.64 


-74.6 


0.082 


0.26 


46 


Continuum 




225.67836 


3.01 x 2.54 


-80.0 


5.330 




0.39 



a Position angle of the long axis of the beam, measured east (counterclockwise) from north. 
" Width of each channel in km s -1 . 

c For lines, the mean of the lcr rms of the spectrum at each spatial position with the spectral resolution given 
in the previous column. For the continuum, the lcr rms of the continuum intensity. 



most likely driving source. Furthermore, as we will show 
below in M5. 11 the position and morphology of CPOC 19 
indicate it is unlikely to be associated with Per-Bolo 58. 

Enoch et al. (2010) obtained deep 70 fxm observations 
of this object as part of a program attempting to detect 
the FHSC by targeting cores classified as starless but 
meeting various criteria suggesting they may already be 
collapsing. Per-Bolo 58 was detected with a low 70 /im 
flux density (65 mjy) consistent with theoretical pre- 
dictions for FHSCs. They also noted an associated 4a 
24 /im detection in existing, shallower Spitzer c2d data 
first published by J0rgensen et al. (2006) 8 above the flux 
level predicted for a FHSC surrounded by a spherically 
symmetric envelope. From continuum radiative transfer 
models, Enoch et al. determined this source has an inter- 
nal luminosity of L, n < ~ 0.012 L©. They showed that the 
observed SED could be fit by either a very low luminos- 
ity protostar or a first hydrostatic core with a spherical 
and/or outflow cavity to reduce the 24 /im opacity, and 
were unable to rule out the possibility that Per-Bolo 58 
is a low luminosity protostar based on the SED alone. 

3. OBSERVATIONS 

One track of observations of Per-Bolo 58 were obtained 
with the SMA on 2010 November 16 in the compact con- 
figuration with 8 antennas, providing projected baselines 
ranging from 5 — 76 m. The observations were obtained 
with the 230 GHz receiver and included 4 GHz of band- 
width per sideband, with 12 GHz spacing between the 
centers of the two sidebands. The correlator was con- 
figured such that the lower sideband (LSB) covered ap- 
proximately 216.8 — 220.8 GHz while the upper sideband 
(USB) covered aproximately 228.8 - 232.8 GHz, provid- 
ing simultaneous observations of the 12 CO, 13 CO, C ls O 
J = 2 - 1, N 2 D+ J = 3 - 2, and SiO J = 5 - 4 lines. The 
lines were observed with either 256 ( 13 CO, SiO) or 512 
( 12 CO, C 18 0, N 2 D+) channels in the 104 MHz bands, 
providing channel separations of 0.53 and 0.26 km s _1 , 
respectively. The remaining bands were used to measure 
the 1.3 mm continuum with a total bandwidth of 5.33 
GHz. 

The weather conditions were extremely good for 230 
GHz observations. Throughout the observations the 
zenith opacity at 225 GHz varied between ~ 0.05 — 0.1, 

8 This 4cr 24 (im detection was not reported in the list of 
candidate young stellar objects (YSOs) reported by j0rgensen et 
al. (2006), but was present in the full source catalog assembled by 
the c2d team (Evans et al. 2007). 



and the system temperature was typically ~ 120K, rang- 
ing from 80 — 240 K depending on elevation. Regular 
observations of the quasar 3C84 were interspersed with 
those of Per-Bolo 58 for gain calibration. Uranus and 
3C84 were used for passband calibration, and Uranus was 
used for absolute flux calibration. We estimate a 15% 
uncertainty in the absolute flux calibration by compar- 
ing the measured fluxes of the calibrators from our cali- 
brated data with those in the SMA calibrator database 9 
for the same observation date. The data were inspected, 
flagged, and calibrated using the MIR software package 10 
and imaged, cleaned, and restored using the Multichan- 
nel Image Reconstruction, Image Analysis, and Display 
(MIRIAD) software package configured for the SMA 11 . 
Imaging was performed with a robust uv weighting pa- 
rameter of +1, found to provide the best compromise 
between resolution and sensitivity for these data. The 
final maps were re-gridded onto 0.5" pixels. 

Only the continuum and CO (and its isotopologues) 
are detected and discussed in this paper. Table Q] lists, 
for each of these observations, the frequency of obser- 
vation, the synthesized beam size and orientation, the 
total bandwidth, the channel separation (for the lines), 
and the measured la rms. For the continuum, the lcr 
rms is determined by calculating the standard deviation 
of all off-source pixels. For the spectra, the la rms is 
determined by calculating, for each pixel, the standard 
deviation of the intensity in each spectral channel outside 
of the velocity range — 20 km s" 1 , and then calculating 
the mean standard deviation over all pixels. 

4. RESULTS 

4.1. Continuum 

Figure [T] displays, in inverted grayscale, the Spitzer c2d 
24 /im image of Per-Bolo 58, showing the weak detection 
at this wavelength. Overlaid are contours showing the 
1.3 mm SMA continuum intensity. Per-Bolo 58 is clearly 
detected in the continuum. Using the MIRIAD task imfit 
we fit an elliptical Gaussian to the detection and report 
a peak position of R.A. = 03 29 25.46, Deck = +31 28 
15.0, as listed in Tabled The uncertainties in this po- 
sition from the fit are less than 0.2" in both Right As- 
cension and Declination. As evident from Figure [TJ this 
position shows excellent agreement with the position of 
the infrared source detected by Enoch et al. (2010). It 

9 Available at http://smal.sma.hawaii.edu/callist/callist.html 

10 Available at https:/ /www. cfa. harvard. edu/^cqi/mircook. html 

11 Available at http://www.cfa.harvard.edu/sma/miriad/ 
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Fig. 1. — Spitzer c2d 24 fim image of Per-Bolo 58. The grayscale 
is inverted and displayed in a linear stretch with the minimum and 
maximum intensities set to 36.5 (white) and 37.1 (black) MJy sr — 1 , 
respectively. The green contours show the SMA 1.3 mm contin- 
uum intensity. The solid contours show positive intensity contours 
starting at 3a and increasing by la and the dotted contours show 
negative intensity contours starting at —3a and decreasing by la, 
where the la rms in the continuum image is 0.39 mjy beam - 1 . The 
two 3cr peaks southeast and northwest of Per-Bolo 58 are sidelobes 
incompletely removed by the cleaning process. Only one — 3a peak 
is present in the northwest corner of the image. The synthesized 
beam size and shape of the SMA continuum observations is shown 
by the white filled oval in the lower right. 

also agrees to within 1.5" with the position reported by 
Schnee et al. (2010) from their 3 mm CARMA detec- 
tion. Our measured peak position is offset by 3.3" to 
the southwest of the Bolocam 1.1 mm single-dish contin- 
uum position reported by Enoch et al. (2010). As the 
Bolocam data includes an absolute pointing uncertainty 
of 7", we do not consider the offset to be significant. 

Other properties derived from the elliptical Gaussian 
fit to the detected source, including the peak flux den- 
sity, the total flux density, the deconvolved source size, 
and the deconvolved source position angle, are also re- 
ported in Table [21 The total flux density and effective 
radius 12 of the continuum source are 26.6 ± 4.0 mjy and 
2.3 ± 0.2" (575 ± 125 AU at the assumed distance of 
250 ± 50 pc), respectively. For comparison, Enoch et 
al. (2006) measured a total flux density and effective ra- 
dius of 330 ± 30 mjy and 21 ± 1" (5250 ± 1080 AU at the 
assumed distance of 250 ± 50 pc) , respectively, in their 
single-dish Bolocam 1.1 mm continuum data. Scaling 
this 1.1 mm flux density to 1.3 mm assuming optically 
thin emission with an opacity power- law index of = 1.5 
gives an expected 1.3 mm flux density of 171 ± 16 mjy 13 . 
This factor of 6 discrepancy in flux density and factor 
of 9 difference in effective radius between the SMA and 
Bolocam observations is easily explained by the fact that 

12 The effective radius is defined as the geometric mean of the 
semimajor and semiminor axes 

13 The uncertainty in the expected 1.3 mm flux density includes 
only the uncertainty in the Bolocam flux density and does not 
include any component from the uncertain dust opacity power-law 
assumption. Power-law indices of 1.0 and 2.0 would give expected 
1.3 mm flux densities of 187 ± 17 and 155 ± 14, respectively. 
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uv distance (kA.) 

Fig. 2. — SMA 1.3 mm continuum visibility amplitude versus uv- 
distance for Per-Bolo 58. The amplitudes are calculated with the 
MIRIAD task uvamp, which bins the data in annuli according to 
uv distance and calculates the vector averaged amplitude in each 
annulus. The annuli are spaced by 5 kA, and the error bars give 
the standard deviation in each annulus. The dotted line shows 
the expectation value for zero signal. The solid line shows the 
best-fit Gaussian function, while the dashed line shows the best- 
fit Gaussian+constant function. In both cases the Gaussian fits 
were calculated using the IDL procedure gaussfit (see text for more 
details). 

less than 3% of the total uv pointings from our compact 
configuration SMA observations are located at projected 
baselines < 10 kA (corresponding to angular scales > 
21"), thus these data filter out the majority of the emis- 
sion from the large-scale extended core and are instead 
sensitive only to the emission from the inner, compact 
envelope and/or circumstellar disk. 

Figure [2] shows the SMA 1.3 mm continuum visibility 
amplitudes for Per-Bolo 58. There may be evidence of a 
compact, unresolved component seen in baselines longer 
than ~ 30 kA. To investigate this possibility, we used 
the IDL procedure gaussfit to fit both a Gaussian and a 
Gaussian plus constant offset to the data assuming Gaus- 
sian errors. As shown in Figure [2 these fits indicate that 
a Gaussian plus constant offset (reduced x 2 = 0.9, with 
a best-fit offset of 4.7 mjy) gives a significantly better 
fit than a Gaussian alone (reduced x 2 — 4.3). We note 
here that if we implement gaussfit with Poisson (y/~N) er- 
rors instead of Gaussian errors, appropriate for the low 
signal-to-noise regime, a Gaussian plus constant offset 
(reduced x 2 — 0.5, with a best-fit offset of 3.3 mjy) still 
gives a better fit by eye than a Gaussian alone (reduced 
X 2 = 1.5), although the quantitative difference between 
the quality of the two fits is small. 

Based on these results we conclude that the presence of 
an unresolved component is possible but only marginally 
supported by the data, especially considering that Gaus- 
sians are likely poor representations of the true emission 
profiles. If real, this unresolved component could repre- 
sent a circumstellar disk or the flattened, disk-like struc- 
ture of a rapidly rotating FHSC (e.g., Saigo & Tomisaka 
2006; Saigo, Tomisaka, & Matsumoto 2008). However, 
at the distance to Perseus, projected baselines ranging 
from 30 — 50 kA correspond to spatial scales ranging from 
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TABLE 2 

Elliptical Gaussian Fit to the 1.3 mm Continuum Detection 



Peak R.A. 


Peak Dccl. 


Peak Flux Density 


Total Flux Density 


Source Size a 


Source P.A. a 


M 


n 


(J2000) 


(J2000) 


(mjy beam -1 ) 


(mjy) 


(") 


(degrees) 


(M Q ) 


(cm -3 ) 


03:29:25.46 


+31:28:15.0 


6.9 ± 1.1 


26.6 ± 4.0 


5.3 ± 0.4 x 4.0 ± 0.3 


19.0 ± 16.2 


0.11 ± 0.05 


2.0 ± 1.6 X 10 7 



a Deconvolved with the beam. 

1000 — 1700 AU, thus the unresolved component could 
also simply correspond to the dense inner envelope. Fu- 
ture observations in more extended configurations pro- 
viding longer baselines must be combined with radia- 
tive transfer models of a star+disk+envelope system fit 
to both the visibility amplitudes and the observed SED 
(e.g., J0rgensen et al. 2005; Enoch et al. 2009b) in order 
to fully separate inner envelope and disk components and 
confirm the presence of this unresolved component. 

The last two columns of Table list the mass and den- 
sity derived from the SMA continuum detection. The 
mass is calculated as 



M = 100 



d 2 S„ 



B u (T d )k v 



(1) 



where S v is the total flux density, B V (T D ) is the Planck 
function at the isothermal dust temperature Td, k v is 
the dust opacity, d — 250 ± 50 pc, and the factor of 
100 is the assumed gas-to-dust ratio. We adopt the dust 
opacities of Ossenkopf & Henning (1994) appropriate for 
thin ice mantles after 10 5 yr of coagulation at a gas den- 
sity of 10 6 cm" 3 (OH5 dust), giving k v = 0.864 cm 2 
gm" 1 at the frequency of the continuum observations. 
For comparison, Enoch et al. (2006; 2010) and Schnee 
et al. (2010) also adopt OH5 dust opacities at the fre- 
quencies of their observations to calculate masses, while 
Hatchell et al. (2007a) adopt a dust opacity 1.5 times 
higher than that of OH5 dust at the frequency of their 
observations. With the total flux density listed in Table 
[5] and an assumed Trj = 10 K, Equation Q] gives a total 
mass of 0.11 ± 0.05 M Q . Note that the uncertainty of 
0.05 M Q only includes the statistical uncertainty in the 
flux density and distance. The true uncertainty is likely 
dominated by the dust temperature and opacity assump- 
tions. For example, the dust opacity at 1.3 mm can vary 
by factors of ~ 2 — 4 depending on which dust opacity 
model is adopted (e.g., Shirley et al. 2005; 2011), di- 
rectly leading to factors of ~ 2 — 4 variation in the mass. 
Furthermore, the exact dust temperature is uncertain. 
Rosolowsky et al. (2008) identified two velocity compo- 
nents in their NH3 observations (see further discussion of 
the two components in ^4. 21 below) with derived kinetic 
temperatures of 10.3 and 10.5 K, respectively, in good 
agreement with our assumption of 10 K. An assumed 
temperature of 20 K would change the derived mass to 
0.04 ± 0.02 Mq. 

The mean number density, n, is calculated assuming 
spherical symmetry as 



M 



(2) 



where M is the mass, r e ff is the effective radius, mn is 
the hydrogen mass, and [i is the mean molecular weight 



per free particle. With M and r e // as given above and 
[i = 2.37 for gas that is 71% by mass hydrogen, 27% 
helium, and 2% metals (Kauffmann et al. 2008), we cal- 
culate n = 2.0 ± 1.6 x 10 7 cm" 3 , where the uncertainty 
includes the statistical uncertainties in the mass and ef- 
fective radius as calculated above. 

4.2. CO and Isotopologues 

Figure H shows the 12 CO, 13 CO, and C 18 J = 2 - 1 
spectra of Per-Bolo 58 averaged over the central 4" rela- 
tive to the position of the continuum source. The C 18 
line is clearly detected as a narrow line centered at ~ 7.3 
km s . Gaussian fits to the individual 0.5" pixels within 
this central region using the IDL procedure gaussfit yield 
central velocities between 7.2 and 7.5 km s _1 , with both 
a mean and median of 7.3 km s . This dispersion in ve- 
locity of 0.3 km s" 1 is approximately within one spectral 
resolution element (0.26 km s _1 ; TableQ]), thus we adopt 
7.3 km s" 1 as the core rest velocity of Per-Bolo 58. For 
comparison, Hatchell & Dunham (2009) derived a core 
velocity of 7.9 km s _1 based on single-dish FCRAO C 18 
J = 1 — data with a spatial and spectral resolution of 
46" and 0.25 km s _1 , respectively, and Kirk, Johnstone, 

6 Tafalla (2007) found a velocity of 7.53 (7.54) km s" 1 
based on single-dish IRAM 30mC 18 OJ = 2- l (N 2 H+ 
J = 1 — 0) data with a spatial and spectral resolution of 
11" (25") and 0.05 km s" 1 (0.05 km s" 1 ), respectively. 
The range of values between 7.3 — 7.9 km s" 1 may be ex- 
plained by the fact that Rosolowsky et al. (2009) argued 
that two velocity components, one at 7.4 and one at 7.77 
km s" 1 , are necessary to fit their NH3 observations. Fur- 
ther detections of dense gas tracers, particularly at the 
high spatial resolution provided by interferometers such 
as the SMA, are needed to further investigate the true 
core velocity. 

The Gaussian fits to the individual pixels within 
the central 4" give full-width half-maximum (FWHM) 
linewidths varying between 0.3 and 1.4 km s" 1 , with 
a mean value of 0.7 km s . The lines are generally 
spread over only 2 — 4 velocity channels and are thus only 
marginally resolved. Purely thermal gas at 10 K would 
produce FWHM linewidths of 0.12 km s _1 and lines con- 
fined to a single channel given the channel width of ~ 0.3 
km s . While nonthermal components thus appear to 
be present, higher density tracers are required to rule 
out contamination by emission from outflowing gas (see 
below). No significant (> 5cr) C 18 emission is detected 
outside of this central region. 

The 13 CO line is also clearly detected, centered at ~ 

7 km s _1 . Gaussian fits to the individual pixels within 
this central region, again using the IDL procedure gauss- 
fit, yield central velocities between 6.8 and 7.3 km s" 1 , 
with both a mean and median of 6.9 km s _1 . We con- 
sider the C 18 line a better measure of the true core rest 
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Fig. 3.— 12 CO, u CO, and C ls O J = 2 — 1 spectra of Per- 
Bolo 58 averaged over the central 4" relative to the position of the 
continuum source. The 13 CO and C ls O spectra are shifted down 
by 1 and 2 K, respectively, for display purposes. 

velocity than the 13 CO line since the latter is more likely 
to include emission from outflowing gas and also more 
likely to be spatially extended and thus partially filtered 
out. We detect 13 CO emission near the core velocity at 
levels greater than 5a over the full extent of our map, 
but the emission is sporadic and not clearly tracing ei- 
ther the central core or outflow axis (see below), and no 
significant emission is detected at \v — v corc \ > 2 km s . 
This emission likely represents the most compact com- 
ponents of extended 13 CO cloud emission, most of which 
is resolved out. Shorter baseline data are necessary to 
fully reconstruct a map. 

Finally, emission from 12 CO is also clearly detected. 
The spectra show two peaks blueshifted and redshifted 
relative to the core velocity and no emission at the core 
velocity, classic signatures of an outflow observed with an 
interferometer (e.g., Bourke et al. 2005) with the emis- 
sion from the core itself resolved out. Figure 2] displays 
channel maps of the 12 CO data at a velocity resolution of 
2 km s . Blueshifted emission between ~ 4 — 7 km s _1 
extends to the east (including in the central velocity 
channel due to the 2 km s _1 channel widths) whereas 
redshifted emission between ~ 10—13 km s _1 extends to 
the west. With our assumed core velocity of 7.3 km s -1 , 
the blueshifted emission appears to be moving slower 
than the redshifted emission. We will return to this point 
in £15.21 No outflow emission is seen beyond |u — u CO re| ~ 7 
km s _1 . 

5. MOLECULAR OUTFLOW 

5.1. Morphology 

Figure [5] shows integrated blueshifted and redshifted 
emission contours overlaid on the 1.3 mm continuum 
image. The east-west outflow is clearly seen, with the 
blueshifted lobe extending almost due east, the red- 
shifted lobe extending almost due west, and very little 
overlap between the lobes. From this observed morphol- 
ogy, it is clear that the candidate outflow feature CPOC 
19 reported by Arce et al. (2010), with redshifted emis- 
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Fig. 4.— SMA 12 CO J = 2 - 1 channel maps of Per-Bolo 58 with 
a velocity resolution of 2 km s — , with the velocities marked in 
the top left of each panel. The gray scale shows intensity in units 
of Jy beam -1 with a linear scaling ranging from 0.0 Jy beam -1 
(white) to 0.50 Jy beam -1 (black). Overplotted with solid lines are 
positive intensity contours starting at 5cr and increasing by 5<r and 
with dotted lines are negative intensity contours starting at — 5<r 
and decreasing by 5<r, where the la rms in the 2 km s — 1 channels 
is 0.022 Jy beam -1 . The synthesized beam size and shape of the 
12 CO J = 2 — 1 observations is shown by the black filled oval in the 
lower right of each panel. The central (0,0) position is that reported 
by Enoch et al. (2010) (the phase center of the observations). 

sion extending to the northeast of Per-Bolo 58, is not re- 
lated to Per-Bolo 58 and is likely a high-velocity feature 
associated with another source in NGC 1333. Enoch et 
al. (2010) noted the presence of faint nebulosity detected 
in the Spitzer Infrared Array Camera (IRAC; Fazio et 
al. 2004) band 2 (4.5 /im) image and suggested this 
emission may be tracing a weak outflow since this band 
overlaps with emission from molecular hydrogen heated 
by outflow shocks (e.g., Noriega- Crespo et al. 2004; De 
Buizer & Vacca 2010). This nebulosity extends a few 
arcseconds to the east of Per-Bolo 58 and is aligned with 
the blue outflow lobe, thus we confirm that it is indeed 
likely associated with the outflow. 

We measure a position angle of 85° (measured east 
from north) for the axis connecting both lobes and pass- 
ing through the position of the infrared source. The 
outflow shows a narrow, jet-like morphology, with open- 
ing semi-angles of ~ 8° for both the blue and red lobes. 
These opening semi-angles are best treated as upper lim- 
its since the lobes are only marginally resolved in the 
direction perpendicular to the flow (total extents of ~ 
7", compared to the 2.9 x 2.5" beam). The blue lobe 
extends 27" (6750 ± 1350 AU at 250 ± 50 pc) from the 
infrared source whereas the red lobe only extends 22" 
(5500 ± 1100 AU at 250 ± 50 pc); we list these values 
in the Table [3] along with other quantities described in 
q!5.2I below. However, given that the infrared source is 
offset by 3.3" from the phase center (i j4.1|) , both lobes 
extend nearly to the half-power radius of the primary 
beam (HWHM - 27" at 230 GHz). The lobe sizes re- 
ported here are thus lower limits to the true sizes. 

Given the morphology of this outflow, in particular 
the two distinct, bipolar lobes showing either entirely 
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Fig. 5. — SMA 1.3 mm continuum image of Per-Bolo 58. The 
grayscale is inverted and displayed in a linear stretch ranging from - 
2.73 mjy beam -1 (—7a, white) to 4.68 mjy beam -1 (12cr, black). 
The green contours show the continuum intensity starting at 3cr 
and increasing by 7cr, where the 1<t rms in the continuum image 
is 0.39 mjy beam -1 . The blue contours show blueshifted 12 CO 
J = 2 — 1 emission integrated from 0.3 to 7.3 km s , while the 
red contours show redshifted emission integrated from 7.3 to 14.3 
km s — 1 . The solid blue and red contours start at 3<r and increase 
by lcr whereas the dotted blue and red contours start at —3a- and 
decrease by lcr, where the lcr rms in both the integrated red and 
blue maps is 0.16 Jy beam -1 km s — 1 . The black circle shows the 
primary beam of the SMA at 230 GHz. The synthesized beam size 
and shape of the 12 CO observations is shown by the black filled 
oval in the lower right. The central (0,0) position is that reported 
by Enoch et al. (2010) (the phase center of the observations). 

blueshifted or redshifted emission, the geometry of the 
Per-Bolo 58 system is best described by case 2 of Cabrit 
& Bertout (1986): the system inclination 14 is greater 
than or equal to the opening semi-angle of the outflow, 
and the sum of the inclination and outflow semi-opening 
angle is less than or equal to 90°. However, with such 
a jet-like morphology with a very small opening semi- 
angle of ~ 8°, we are only able to rule out extreme edge- 
on and face-on inclinations. In general we are unable to 
set strong constraints on the inclination, a point we will 
return to in SjB] below. 

5.2. Kinematics 

We calculate, separately for the blue and red lobes 
and for the combined outflow, the mass of the outflow, 
Mfi ow , and its kinematic and dynamic properties (kinetic 
energy [Efi ow ], momentum [Pfi ow ], characteristic veloc- 
ity [v = Pfiow/Mfi ow ], dynamical time [r d = R lobe /v], 
mechanical luminosity [Lfi ow — Efi ow /Td], and force 
[Ffiow — Pfiow/Td}) in the standard manner (e.g., Cabrit 
& Bertout 1986, 1990; see also Dunham et al. 2010). 
The results are listed in Table [3] using a map corrected 
for primary beam attenuation. The mass (and all other 
properties that depend on mass) are calculated assum- 
ing LTE and do not include corrections for inclination or 
optical depth. Thus, combined with the facts that the 
outflow likely extends beyond the primary beam ( ^5.ip 
and that these observations will filter out any extended 

14 An inclination of 0° corresponds to a face-on system, whereas 
an inclination of 90° corresponds to an edge-on system. 



component to the outflow, our reported values are lower 
limits only. Because we only calculate lower limits, we 
assume the gas temperature that minimizes the LTE cal- 
culation of mass, 17.6 K (see Dunham et al. 2010 for de- 
tails). The actual temperatures of outflowing gas may 
be significantly higher (50 - 200 K; Hatchell et al. 1999; 
van Kempen et al. 2009a, 2009b), although this remains 
a significant unknown, especially for very low luminosity 
sources like Per-Bolo 58. The true values of the quantities 
listed in Table |3] may be higher by more than an order 
of magnitude depending on the true values of the out- 
flowing gas temperature and optical depth in the 12 CO 
J = 2 — 1 line and the fraction of total outflow emission 
missing from our map (see, e.g., Bourke et al. 2005, Chen 
et al. 2010, Arce et al. 2010, Curtis et al. 2010, and Dun- 
ham et al. 2010 for a few of many recent discussions on 
calculating outflow properties). 

While the morphology of the blue and red lobes ar- 
gues strongly for the presence of an outflow, with such 
a low characteristic velocity (2.9 km s _1 ) it is important 
to rule out bound motion as being responsible for the 
kinematic signatures interpreted as an outflow. For gas 
at 2.9 km s -1 to be bound at 6200 AU (the average ob- 
served size of the two lobes) requires a mass of at least 
29 M Q , a factor of ~ 15 — 30 times higher than the total 
core mass. We thus definitively eliminate bound motion 
as being responsible for the observed blue and red lobes. 

The mass of outflowing gas in the blue and red lobes 
are comparable (~lx 10~ 4 M ), but the characteristic 
velocity of the red lobe is faster by 1.1 km s _1 (46%), 
leading to larger values of all other quantities except for 
the dynamical time, which is inversely proportional to 
velocity and thus smaller. In particular, the force re- 
quired to drive the red lobe is a factor of 2 higher than 
that required to drive the blue lobe. Such a force mis- 
match could, in principle, accelerate the central object 
out of the dense core from which it is forming. In reality, 
however, this acceleration is very small; a simple calcu- 
lation assuming a central mass of 0.01 M Q 15 shows that 
it would take 0.54 Myr for the central object to move 
2.3" (575 AU at 250 pc, the effective radius of the SMA 
continuum source). This time is comparable to the to- 
tal duration of the embedded phase (0.44 — 0.54 Myr; 
Evans et al. 2009). Furthermore, if we shift the central 
velocity of the core from our assumed value of 7.3 km s -1 
to 7.9 km s _1 , as given by Hatchell & Dunham (2009), 
the forces required to drive the two lobes agree to within 
16% and the mismatch is negligible. Detections of higher 
density tracers than C 18 0, used by both Hatchell & Dun- 
ham (2009) and by us to determine the core velocity, are 
required to better pin down the true velocity of Per-Bolo 
58. 

6. DISCUSSION 

Between Per-Bolo 58 (Enoch et al. 2010; this work), 
L1448 IRS2E (Chen et al. 2010), and L1451-mm (Pineda 
et al. 2011), three candidate FHSCs have been identi- 
fied in Perseus. Enoch et al. (2009a) identified 66 em- 
bedded protostars in Perseus based on the combination 
of complete surveys of the cloud with the Spitzer c2d 

15 Approximately the amount of mass that would accrete in the 
dynamical time of the outflow (~ 10 4 yr) at a standard rate of 
~ 10 -6 M yr -1 . 



TABLE 3 
Outflow Properties 



Quantity 





Unit 


Blue Lobe 


Red Lobe 


Combined 


AU 




6.8 X 10 3 


5.5 X 10 3 


6.2 X 10 3 


M 




1.1 x 1CT 4 


0.9 x 10- 4 


2.0 x 10~ 4 


M 


km s~ 1 


2.5 x 10~ 4 


3.1 x 10~ 4 


5.6 x 10~ 4 


ergs 




2.3 x 10 40 


2.8 x 10 40 


5.1 x 10 40 


L© 




1.2 x 10~ 5 


2.6 x 10~ 5 


4.2 x lO" 5 


M 


km s — 1 yr 


1.6 x 10 -8 


3.5 x 1(T 8 


5.6 x 10~ 8 


km ; 


,— 1 


2.4 


3.5 


2.9 






1.6 x 10 4 


0.9 x 10 4 


1.0 x 10 4 



Lobe Size (Ri obe ) 
Outflow Mass (M flow ) 
Outflow Momentum (Pfi ow ) 
Outflow Kinetic Energy (Efi on! ) 
Outflow Luminosity (Z//; ou ,) 
Outflow Force {Ffi om ) 
Characteristic Velocity (vfi on; ) 
Dynamical Time (t<j) 



Legacy project (J0rgensen et al. 2006; Rebull et al. 2006) 
and Bolocam 1.1 mm continuum observations (Enoch et 
al. 2006). Assuming the total duration of the embedded 
stage is 0.54 Myr (Evans et al. 2009) and the lifetime of 
the FHSC stage is between 5 x 10 2 — 5 x 10 4 yr (Boss 
& Yorke 1995; Omukai 2007; Saigo, Tomisaka, & Mat- 
sumoto 2008; Tomida et al. 2010), there should be < 1-6 
FHSCs in Perseus. The current number of candidates 
(three) is near the maximum of this range, suggesting 
a relatively long duration of the FHSC stage (at least 
2.5 x 10 4 yr) if all three are in fact first cores. 

With a characteristic velocity of 2.9 km s~ x and a max- 
imum velocity of ~ 7 km s _1 , the Per-Bolo 58 outflow 
is relatively slow compared to other collimated outflows 
from young, embedded objects (e.g., Gueth & Guilloteau 
1999; Lee et al. 2000; Arce & Sargent 2006). Such slow 
velocities agree with the prediction by Machida, Inut- 
suka, & Matsumoto (2008) that the FHSC drives a slow 
(< 5 km s _1 ) outflow. The outflow also features a colli- 
mated, jet-like appearance, with an opening semi-angle 
of ~ 8° and a lobe extent-to-width ratio (E/W) > 3 - 4. 
Machida, Inutsuka, & Matsumoto (2008) predict that 
the outflows driven by first cores have very wide open- 
ing angles and constant E/W ^ 2.2 — 2.5, whereas those 
driven by second cores (protostars) are more collimated 
with E/W starting around 5 and increasing with time. 
The outflow driven by Per-Bolo 58 appears to feature a 
morphology between these two cases but more consistent 
with a protostellar rather than FHSC outflow, although, 
as noted in ^4.11 less than 3% of the total uv point- 
ings from our compact configuration SMA observations 
are located at projected baselines < 10 kA (correspond- 
ing to angular scales > 21"). Future observations pro- 
viding shorter baselines and/or deeper single-dish data 
than published by Hatchell & Dunham (2009), Curtis et 
al. (2010), and Arce et al. (2010) should be pursued to 
determine whether or not a wider component to the out- 
flow more extended than ~ 21" exists and is resolved out 
in the current SMA data. 

Table [3] lists the outflow properties of several other 
embedded sources with similarly low luminosities. L1448 
IRS2E is a dense core driving a molecular outflow with no 
associated Spitzer infrared source and has been identified 
as a candidate FHSC by Chen et al. (2010). Similarly, the 
dense cores CB 17 MMS (X. Chen et al. 2011, in prepara- 
tion) and L1451-mm (Pineda et al. 2011) are associated 
with outflow detections (tentative in the case of CB 17 
MMS), lack associated Spitzer infrared sources, and have 
been identified as candidate FHSCs. L1014-IRS is a very 



low luminosity (~ 0.09 L Q ; Young et al. 2004) Class 
protostar driving a compact outflow detected only in in- 
terferometer observations (Bourkc et al. 2005). Similarly, 
L1148-IRS is a borderline Class 0/1 protostar with Li nt 
~ 0.10 L Q also driving a compact outflow detected only 
in interferometer observations (Kauffmann et al. 2011). 
Finally, L673-7-IRS is a very low luminosity (~ 0.04 L©) 
Class protostar that is driving a large-scale outflow 
(Dunham et al. 2010). Figure [5] plots these various out- 
flow properties versus L int for these six sources plus Per- 
Bolo 58. Significant dispersion of at least one order of 
magnitude (and sometimes much more) is present in all 
properties, and no obvious correlation with luminosity 
is seen. Per-Bolo 58, CB 17 MMS, and L1451-mm fea- 
ture slow outflow velocities consistent with predictions 
by Machida, Inutsuka, & Matsumoto (2008) for FHSC 
outflows, whereas L1448 IRS2E drives an outflow with 
a characteristic velocity (25 km s -1 ) much faster than 
predicted for first cores. However, the protostars L1014- 
IRS, L1148-IRS, and L673-7-IRS also drive slow (< 3.5 
km s _1 ) outflows, and in general there is no clear sepa- 
ration between the properties of the outflows driven by 
the candidate FHSCs and by the very low luminosity 
protostars. 

While these findings argue against a distinct difference 
in evolutionary status between these two groups of ob- 
jects, we caution that the sample size is extremely small, 
three of the four candidate FHSCs lack any mid and 
far-infrared detections and thus only have upper limits 
for their luminosities, and the outflow properties them- 
selves are highly uncertain due to nonuniform assump- 
tions about outflowing gas temperature, source inclina- 
tion, optical depth, and methods of calculating prop- 
erties. This figure must be revisited as more sources 
(both candidate FHSCs and low luminosity protostars) 
with well-constrained inclinations are identified and their 
luminosities and outflow properties are calculated from 
uniform datasets with consistent methods, and as addi- 
tional simulations fully explore predicted differences be- 
tween the kinematics, energetics, and morphologies of 
outflows driven by first cores and protostars. 

A significant unknown in our analysis of the kinematics 
and morphology of the Per-Bolo 58 outflow and compar- 
ison to the predictions of Machida, Inutsuka, & Mat- 
sumoto (2008) is the source inclination. As discussed 
above, we are unable to set strong constraints on the in- 
clination from the outflow morphology. Both relatively 
edge-on and relatively face-on configurations are possi- 
ble, with only extreme edge-on (> 80°) and face-on (< 
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TABLE 4 

Comparison of Outflow Properties 



Quantity 


Per-Bolo 58 


L1448 IRS2E a 


CB 17 MMS b 


L1451-mm c 


L1014-IRS d 


L1148-IRS C 


L673-7-IRS f 


Source Internal Luminosity (£j n t; Lq) 


0.012 


<0.1 


<0.07 


<0.03 


0.09 


0.10 


0.04 


Lobe Size (-Ri f, e ; AU) 


6.2 x 10 3 


9.6 x 10 3 


8.0 x 10 3 


5.5 x 10 2 


5.4 x 10 2 


1.8 x 10 3 


4.5 x 10 4 


Outflow Mass (Mfi ow ; M Q ) 


2.0 x 10" 4 


2.0 x 10" 3 


7.9 x 10~ 4 


1.2 x 10~ 5 


1.4 x 10- s 


8.5 x 10~ 4 


5.0 x 10" 2 


Outflow Momentum (Pf; ou ,; Mq km s — 4 ) 


5.6 x 10~ 4 


5.0 x 10" 2 


1.9 x 10~ 3 


1.7 x 10~ 5 


2.4 x 10~ 5 


8.5 x 10~ 4 


1.3 x 10" 1 


Outflow Kinetic Energy (Efi ow ; ergs) 


5.1 x 10 40 


1.2 x 10 43 


4.3 x 10 40 


3.1 x 10 38 


4.8 x 10 38 


8.0 x 10 38 


6.1 x 10 42 


Outflow Luminosity (Lf; om ; Lq) 


4.2 x 10" 5 


5.0 x 10" 2 


2.5 x 10~ 5 


1.3 x 10~ 6 


3.4 x 10" 6 


8.5 x lO" 6 


8.3 x 10~ 4 


Outflow Force (Fa ow ; Mq km s _1 yr _1 ) 


5.6 x 10" 8 


2.5 x 10" 5 


1.3 x 10" 7 


8.3 x 10" 9 


3.4 x 10" 8 


1.0 x 10" 7 


2.1 x 10" 6 


Characteristic Velocity (v{i ow ; km s — 1 ) 


2.9 


25 


2.4 


1.3 


1.7 


1.0 


3.5 


Dynamical Time (r^; yr) 


1.0 x 10 4 


1.8 x 10 3 


1.4 x 10 4 


2.0 x 10 3 


7.0 x 10 2 


8.5 x 10 3 


6.1 x 10 4 



Cho 



ct al. (2010) 



X. Chen ct al. (2011, in ; 

Pineda ct al. (2011) 
1 Bourkc ct al. (2005) 
! Kauffmann ct al. (2011) 
Dunham ct al. (2010) 
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Fig. 6. — Various outflow properties plotted vs. Li„t for the sources listed in Table|4] Filled black symbols denote the candidate FHSCs 
and filled gray symbols denote the low luminosity protostars. Each source is labeled, and the triangles pointing to the left denote upper 
limits in luminosity for three of the four candidate FHSCs (L1448 IRS2E, CB 17 MMS and L1451-mm). Note that these upper limits are 
actually upper limits to Li, i rather than Li n t and that L{ n t is likely less than L(, D ; , which will include a component from external heating. 
However, without any far-infrared detections, it is not possible to separate internal and external components to L(, ; for these sources. 



10°) inclinations ruled out 16 . If Per-Bolo 58 is viewed 
along a relatively edge-on line-of-sight, the true velocity 
of the outflow would increase beyond that expected for 
FHSCs. For example, an inclination of 80° would in- 
crease the outflow velocity to 16.5 km s . On the other 
hand, if Per-Bolo 58 is viewed along a relatively face-on 
configuration, the true outflow velocity would be similar 
to the characteristic radial velocity of 2.9 km s _1 and 

16 Enoch et al. (2010) also ruled out inclinations > 80° based on 
their radiative transfer models. 



the observed SED would only be consistent with FHSC 
models and not protostellar models (Enoch et al. 2010). 
However, the extent-to- width ratio would increase sub- 
stantially beyond that predicted by Machida, Inutsuka, 
& Matsumoto (2008) for FHSC outflows and the dynam- 
ical time would increase beyond the maximum expected 
lifetime of first cores. For example, at an inclination of 
10° the outflow velocity would remain at 2.9 km s" 1 but 
the lobe size would increase to 3.5 x 10 4 AU, the extent- 
to- width ratio would increase to ~ 20, and the dynamical 
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time would increase to 5.7 x 10 yr. 

Clearly, knowledge of the source inclination would aid 
in determining the evolutionary status of Per-Bolo 58. 
Inclination constraints could be provided by very deep 
near- or mid-infrared images that detect scattered light 
cones from the outflow cavities, combined with radia- 
tive transfer model images at different inclinations (e.g., 
Huard et al. 2006; Terebey et al. 2006). They could also 
be obtained through very high-resolution interferometer 
observations separated by several years that measure the 
proper motion of the outflow (e.g., Lee et al. 2009). 

If Per-Bolo 58 turns out to be a protostar rather than 
a FHSC, our current understanding of the protostcllar 
luminosity distribution would be incomplete. Evans et 
al. (2009) and Enoch et al. (2009) identified 112 embed- 
ded protostars in the Spitzer c2d Legacy dataset and 
showed that their luminosities vary between ~ 0.05 — 50 
L©, with the distribution skewed towards low luminosi- 
ties (59% have L < 1.6 L Q ). In a closely related investi- 
gation, Dunham et al. (2008) presented a detailed study 
of the low end of the protostellar luminosity distribution 
and found that the distribution, which generally rises to 
lower luminosities, flattens out below 0.1 L Q and appears 
to remain flat down to their sensitivity limit of 0.004 L Q 
at 140 pc (or 0.013 L© at 250 pc, the distance to Perseus). 
They noted that the fact that sources were found all the 
way down to their sensitivity limit implied there could 
be a population of objects at even lower luminosities. 
Per-Bolo 58 as a protostar would be a member of such a 
population. Future work must address this issue by fur- 
ther study of the low end of the luminosity distribution of 
embedded sources through a combination of sensitive in- 
terferometer outflow surveys and very deep Herschel and 
James Webb Space Telescope infrared surveys directed 
towards cores currently classified as starless to identify 
other objects similar to Per-Bolo 58, L1448 IRS2E, CB 
17 MMS, and L1451-mm. 

7. SUMMARY 

In this paper we have presented new 230 GHz Submil- 
limeter Array observations of the candidate first hydro- 
static core Per-Bolo 58 and reported the detection of a 
continuum source and a bipolar molecular outflow, both 
centered on the position of the infrared source embedded 
within this core. 

The continuum detection has a total flux density of 
26.6 ± 4.0 mJy, and from this detection we calculate a 
mass of 0.11 ± 0.05 M and a mean number density of 
2.0± 1.6 x 10 7 cm~ 3 . The total mass of the Perseus Bolo 
58 core is approximately 1 M Q , indicating that the SMA 
filters out all but the compact, dense inner region. Fits 
to the visibility amplitudes suggest the presence of an un- 
resolved component with a total flux density of 4.7 mJy, 
but observations in more extended configurations provid- 
ing longer baselines are needed to confirm the presence of 
this unresolved component and determine whether it is 
tracing the dense inner envelope or a circumstcllar disk. 



The outflow is detected in the 12 CO J = 2 — 1 obser- 
vations and extends along a nearly due east- west axis at 
a position angle of 85° (measured east from north) . The 
blue lobe extends to the east and the red lobe extends to 
the west with very little overlap between the lobes. The 
outflow is slow (characteristic velocity of 2.9 km s _1 ), 
shows a jet-like morphology (opening semi-angles ~8° 
for both lobes), and extends to the edges of the pri- 
mary beam. The morphology and kinematics of this out- 
flow are partially consistent with predictions for outflows 
driven by both first hydrostatic cores and protostars, and 
do not fully distinguish between the two possibilities. We 
calculate the kinematic and dynamic properties of the 
outflow in the standard manner and compare them to 
several other protostars and candidate first hydrostatic 
cores with similarly low luminosities, and we show there 
are no clear trends between these properties and source 
luminosity, and no clear separation in properties between 
outflows driven by candidate first hydrostatic cores and 
very low luminosity protostars. 

While Per-Bolo 58 is a very strong candidate to be 
a first core, ultimately we are still unable to determine 
the true evolutionary status of this source; it could be 
either a first hydrostatic core or a very low luminosity 
protostar. We have discussed the evidence in support of 
both possibilities, as well as the implications for our cur- 
rent understanding of protostellar luminosities if it is in 
fact a protostar. We have outlined future observations 
needed to better determine the physical and evolutionary 
properties of this source, including longer baseline con- 
tinuum data and source inclination constraints provided 
by either deep infrared images or very high-resolution 
interferometer observations, and we have also discussed 
the need for further study of the low end of the protostel- 
lar luminosity distribution through both sensitive inter- 
ferometer outflow surveys and very deep Herschel and 
James Webb Space Telescope infrared surveys directed 
towards cores currently classified as starless. 
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